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Abstract 

The driving forces for a generally oriented crack embedded in a Functionally 
Graded strip sandwiched between two half planes are analyzed using singular integral 
equations with Cauchy kernels, and integrated using Lobatto-Chebyshev collocation. 
Mixed-mode Stress Intensity Factors (SIF) and Strain Energy Release Rates (SERR) are 
calculated. The Stress Intensity Factors are compared for accuracy with previously 
published results. Parametric studies are conducted for various non-homogeneity ratios, 
crack lengths, crack orientation and thickness of the strip. It is shown that the SERR is 
more complete and should be used for crack propagation analysis. 


1. Introduction 

One way to reduce the residual stresses in composites is to process fully tailored 
materials and interfacial zones with predetermined continuously varying mechanical 
properties. Such materials are known as Functionally Graded Materials (FGM) (see 
Asish et. ah, 1997 and Holt et. ah, 1993). Some FGM could be described as two-phase 
particulate composites where the volume fractions of its constituents differ continuously 
in the thickness direction (see Niino and Maeda, 1990; Hirano and Yamada, 1988; Hirano 
et. ah, 1988; and Kawasaki and Watanabe R., 1990). This implies that the composition 
profile could be tailored to give desired thermomechanical properties. One of the most 
important of these properties is the minimization of crack propagation. In order to design 
FGM components, then, the driving forces of crack propagation must be fully 
understood. 

The problem under consideration here is that of a generally oriented crack embedded 
in a nonhomogeneous strip sandwiched between two isotropic half planes. A system of 
singular integral equations with Cauchy kernels is used to analyze the driving forces 
(Stress Intensity Factors (SIF) and Strain Energy Release Rates (SERR)) of crack 
propagation. 

The present work is a generalization of a sequence of papers (Delale and Erdogan 
(1983), (1988a), (1988b), Konda and Erdogan (1994), and Chen and Ergodan (1996)) 
concerning driving forces of crack propagation for problems involving various boundary 
conditions and crack geometry. In these papers, an exponential variation in material 
properties within the FGM is assumed, and it is shown that Poisson's ratio has little effect 
on stress intensity factors. Therefore, in the present formulation, the same Poisson ratio 
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is used in all three materials, while the shear modulus has an exponential form. Also, 
these papers considered only horizontally oriented cracks, while the present paper 
addresses a crack with arbitrary orientation angle. 

The solution methodology follows the basic steps presented in the previous papers. 
Specifically, the problem is cast in perturbation form. First, the crack surface tractions 
are computed for an FGM embedded between the two half planes with given far field 
stresses when no crack is present. In the second step, these tractions are used to compute 
the stresses at the crack tips for the perturbation problem. These steps are depicted in 
Figures lb and lc. In order to account for the arbitrary orientation angle, the perturbation 
problem is separated into two parts, depicted in Figures 2b and 2c. The first part includes 
the influence of the interfaces, and the second part examines the crack in an infinite 
FGM. Details are shown in the next section. 


2. Formulation of the problem 

The geometry of the problem is shown in Figure la. The two dissimilar materials, 
which are perfectly bonded to the FGM, are isotropic and homogeneous, the FGM has a 
finite thickness h, and is denoted as Material 2. Material 1 occupies the lower half plane, 
for y < 0, while Material 3 occupies the upper half plane for y > h. 

In global coordinates (x,y), the shear modulus of the FGM is assumed to be as 
follows: 

// 2 (y ) = fi l e ry (1) 

and in local coordinates (xi,yi) as: 

/n 2 (x„y\) = H i e Px ' +s " ( 2 ) 

where 

Y = ~ln(— ) 

h jU, 

5 = ycos(d) (3) 

p = ysin(0) 

Hooke’s law relates strain and stress using two independent material constants: 


e*x =7 -[(k + 1)o„ +(k-3)o ¥V ] 

8p 


£ y> +( K+1 ) a yy] 


£ 


xy 



(4) 
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where the bulk modulus, k, is defined as. 


k = 3-4v for 


plane strain 


or 


3-v 

K = 

1+v 


for plane stress 


(5) 


( 6 ) 


The solution strategy is shown in Figures lb and lc. The governing equation for the 
half plane is 


<9 4 F,(A,y) | ^ d 4 F\(x,y) | <9 4 F,(A,y) _ 
dx 4 “ dx-dy 2 d y 4 

The solution of (7) is found by applying the Fourier Transform: 


F, (a-, y) = J V (a, y)e " a da 


(7) 


( 8 ) 


and solving the characteristic equation 


4 ., . ,d 2 V d 4 V n 

a V-2a" r + r = 0 

dy^ dy 4 

==> m 4 -2a : m : +a 4 =0 
=> m, = m, = loci & m 2 = m 4 = 


Ha 


(9) 


so that. 


F, (x, y) = -^f [(D, (a) + yD, (a))e |a|y + (D,(a) + yD 4 (a))e >|y ] e lxa da (10) 
2tc j 


Due to the condition of irregularity at y<0 (the stress function vanishes as y 
D^oO and D 4 (a) must be zero, therefore 


-oo), 


F, (A, y) = — f [(D, (a) + yD, (a)k 1 " 1 ' \e" a da 
2 K J 


(i i; 


Similarly for Material 3, the application of the condition of irregularity at y>0 yields 
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( 12 ) 


F 3 ( x, v) = — f[(C, (a) + yC 2 (a))e M ' \e" <, da 
In J 


The solution strategy of the perturbation problem is shown in Figure 2. The Airy 
stress function method is adopted in this study, mainly for making use of the technique 
developed by Delale and Erdogan (1988). It is assumed that the Airy stress function for 
the FGM is composed of two functions; one is associated with an infinite plane 
containing the crack on the x r axis, U 2 (Xi,yi) (see Figure 2c), while the second is an 
uncracked strip, F 2 (x,y) (see Figure 2b). 

The governing equation in the global coordinate system is 


V 2 (c xx +O yy ) + y 2 
V 4 F 2 (x,y) + y 2 


k, -3 

+~” 7°v 
K ? + 1 


-2y— (c xx+ c vy ) = 0 
dy 


d 2 F 2 (x,y) K, -3 d 2 F,(x,y) 


dy 2 


K, +1 dx' 


(13) 


-2y— V 2 F,(x,y) = 0 
dy 


The characteristic equation of ( 13) is 


m 4 -2ynr +(y 2 -2a 2 )nr +2ya : m + (a 4 


ay 


k, -3 

K 0 + 1 


■) = o 


(14) 


Four roots of ( 14) and the stress function are obtained in the following forms: 


m, 




, - k 2 

V K, +1 


y y 2 7 3- k, 

m, = — + Ja" + — + iay 

- 2 \ 4 \ k 2 +1 


m, = 


1 

2 



y r y 2 " /3 -k 7 

2 \ 4 V k 2 + 1 


(15) 


F,(x,y) =— f[A,(a)e ra,y +A,(a)e m2y +A 3 (a)e m,y + A 4 (a)e m4> ] e‘ 

2k j 


"da 


(16) 


The governing equation in the local crack coordinate system is 
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( 17 ) 


V 4 U 2 (x,,y l )-25-^— (V 2 U 2 (x 1 ,y,))-2P— — (V 2 U 2 (x l ,y,)) + 


3y. 


dx, 


8p8 d 2 U 2 (x,,y,) | 2 ^ k 2 -33 2 U 2 (x,,y,) | 9 2 U 2 (x,,y,) ) | 


1+k 2 3x ,3y , 


1 + k 2 3yf 


dxf 


52( a 2 U 2 (x,,y,) | k 2 — 3a 2 U 2 (x,,y,) ^ = Q 


8yf 


1 + K, 


dx: 


The characteristic equation of (17) is 


n 4 - 28n 3 + 


P 2 — — - + 8 2 -2a(i(3 + a) 


Kt + 1 


n 2 + a8 


8P ; 


K n + 1 


i + 2a\i 


+ or 


3-k, 


a 2 + 8 : — — — + P(2ia — P) 

K, + 1 


= 0 


(18) 


The four roots of ( 1 8) are 


n, =■ 


n, = 


" 3 = 9 


s+P, 


s-p. 


8 + P, 


s-P, 


| 3- K, 

K, + 1 


l_ 

2 


3 - K, 
K, + I 


j 3-k 2 

K, + 1 


3 - k 2 
K, +1 


8 + P, 


3- k 2 

K, + 1 


( 


+ 4(or +ia 


P-8, 


3-k 2 
'k, +1 


8- Pi 


' 3- K, 

K, + 1 


V 


+ 4(a 2 +ia 


P + 8, 


3-k 2 
K, + 1 


H 

? 


+ - 


8 + P, 


3- k 2 

K, + 1 


V 


( 


+ 4(or +ia 


P-8, 


3- k 2 

K, + 1 


8-P, 


3- k 2 

K-, + 1 


+ 4(or +ia 


P + 8, 


3-k 2 
'k, + 1 


(19) 


By examining the roots carefully, it can be noticed that ni and m are always negative 
as a— >+», while n 3 and n 4 are always positive as a— »±°°. This implies that the stress 
function can be expressed in the following form: 


1 

U, (x. , v, ) = — f[S. (a) e' h> ' + B, (a)e" 2 '' ] e"'"da\ 
" 2 n i, 

£/,(jr 1 , v, ) = — f [5, (a) + B 4 (a)e’ u " ] e‘" a da\ 

In J 


v, > 0 


y, < 0 


( 20 ) 
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Two constants in the system of equation (20) can be determined by application of the 
continuity of stresses at yi=0, as follows: 

3U,(x 1 ,0 + ) _ aU 3 (x,,Q-) 

3y, dy, (21) 

U 2 ( x , ,0 + ) = U 2 (x , ,0” ) 

=> B,(oc) = — — — B,(a)+— — — B.(a) 
n 4 -n ? n 4 -n ? ' 

B 4 (a) = B, (a) + B 7 (a) 

n 4 -n, n 4 -n, 

The general forms of the stress functions used to generate the stress components due 
to each problem have been obtained. Next, formulate the stresses for the infinite plate 
with a crack by differentiating the stress function IN (B ■ and B 2 are still functions of a): 


_ KM, +. 

Ur.V, ) = 


2 K 


j[n;B t e"'" +n;B 2 e* :y '] e ix ' a da 


~ Jcr[ fl/"' 1 + 


da 


cr i i a ■ ’(.q, v, ) = — f[/j*(tv,B, + vv,2?,)e"'- Vl + n] ( h’ 4 B, + w 4 Z?, )c" J ' 1 1 e"' a da (23) 
11 2 n J 


— 1 °° 

a, v " : 'U,, y,~) = — farin',#, + w,B,)er" , - Vl + (vv\B ] + w^B^e"^'] e"' a 
11 2k J 


da 


th 11 1 f ^\ 

’ — — ! — • » 4 1 — — ± _ * U = — : * W — 


tu - n 


-;hs = — — -;w\ = 

- n A - n A 


Hy - n ] 


= 


ru, - « 2 

W 3-04 


The singular integral equations for this class of problems are formulated in terms of 
two auxiliary functions: 


f,(x,) = 2— [u l (x 1 ,0 + )-u l (x l ,CT)] 
dx, 

f,(x l ) = -^-[v l (x r 0 + )-v l (x l ,0")] 
dx, 


(24) 


Notice that the auxiliary functions are valid for any X|, but are nonzero only within the crack (a. b). 
Using Hooke's law for the stresses given by (23), local strains and displacements can 
be calculated. For example, 
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( 25 ) 


e^ix^yl) = -—2- J[(k% + 1 )(nf + n;B.e '"" ) 

8 p 2 2 ni ' - - 

-a : (K 2 -3)(B t e’ hy + B 2 e n - y ' )]e‘ x,a da 

£, U ' Up)) ) = ~ — — — ?[(k? + l)(n;(w i B ] + w 1 B 1 )e"'' 1 + n 4 (vv\B i + w 4 B-,)e n * >l ) 

8 /t-, 2 k ^ ' 

-a : (K--, -3)((u’,5 1 + w 2 B 2 )e"' y ' +(w\B l + w 4 B 2 )e’ uy ' )]e' x ' a da 


Using (25) and first equation of (24) it can be shown that 


1 

8/v^/iU'i) = — [(^2 +l)[(»i 2 -nj w, -n]w 2 )B s 

sir J 


+ (n; — n 2 w 2 - n] w’ 4 )B 2 ] e"' a da 


The Fourier transform of (26) yields 


F t {oc) - /I, , fi, + h {2 B 2 

K ~hl 

/i M =-2 (// -«,)(«, -« 4 ) 

8 /t, 

+ 1 

=“f {n, -n^){n,-n 4 ) 

F 1 (a) = J/ 1 (r) i ' ( ^' a Vr 


In order to determine f 2 (xi), it is necessary to find v(X|,yi) by integrating the normal 
strains in the y direction: 


■u- , y ; ) = - 3)-l- ? B x e'"" + B,e'"" )e'" a da 

8/t, " 2K^n t -o n 2 -o 

- (k 2 + 1)2- la 2 (— — Bye™ + — ’ — )^' V| '^a] 

/ 7 T » I? — A u — A 




e”> yi +—^B A e >ux ' )e l " a da 


8 // j “ 2;r n 3 - <5 n 4 - <5 


-(*, +D— [cr( -S 4 ^' hV| + -B+e^'y^da] 

2 k ^ « 3 - o n 4 - o 
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The constant of integration can be set to zero due to the fact that the plate is fixed at 
the origin. Differentiating (28) and substituting into second equation of (24), then taking 
the Fourier transform yields 


F 2 {a) = h 2] B ] +h 22 B 2 

h, t = ia o - {a 2 (l + K,)-S 2 (K, -3))[ 
8/i, 

h „ = lCX - {a 2 { [ + Ki)-S 2 (k j -3))[ 

8/r, 


(», -^Xn, -h 4 ) 

(5 -n^iS -n 2 )(S -n 4 ) 
(n 2 -n, )(/?-, -n A ) 
(S-n 2 )(S-n 3 )(S-n 4 ) 


] 

] 


b 


F 2 (a) = jf 2 (t)e {p -‘ a), dt 

a 


(29) 


Equations (27) and (29) can then be solved for Bi and Bi: 

F,(a)h 22 -F 2 (a)h l2 

D 1 

h 1 1 h 22 — h 21 h 12 
-F,(a)h 21 + F 2 (cc)h n 
h 1 1 h h^|hjo 


(30) 


The stresses at any given point in the cracked FGM strip can be expressed by the 
sum of stresses obtained from the U 2 and F 2 Airy stress functions, namely: 


(2i 


/ v _ (((' 1 ) / 

U,y,) = ( 7 „ - (jfi.v 




for (/, j = x, , v, ) 


(31) 


These are expressed in (x,y) or (x t ,yi) coordinates using the regular stress transformation: 


f 2 

<*.u 


m 2 n 2 — 2mn 




► = 

n 2 m~ 2mn 


CT v,v, ► 






T ". 


mn — mn in —n~ 


. T ,,v, . 


m = cos(0);aj = sin(0) 


(32) 


thus the stresses for the FGM are. 
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G <;’(x,y) =— f[mfA,e m ' y +m;A,e m ^ y +m 1 2 A,e n ’’ y +m^A 4 e n,J> ]e i 
2k J 


x da 


+ cos(0) 2 y- J[n 2 B,e n ' yi +n 2 B 2 e" :y ' ]e 1X| “da 


f <2> 

yy 


-sin(6) 2 — for[B,e n ' y ' + B.e"^ ]e' m da 
2k _{ 

1 °° 

+ 2sin(0)cos(0) — | ia[n l B,e n ' y ' +n,B 1 e" !!l ]e"‘ l "da 
2k _ j x 

;»(x,y ) = -— [a 2 [A,e m,y + A,e m:y + A 2 e m,y + A 4 e n,jy ]e ixa da 
2k 

+ sin(0) 2 — T[n 2 B l e n ' y ' +n;B^e" :>l ] e lx ' a da 
2k J 


— cos(0) 2 — — f a 2 [B,e n ' yi + B,e n;>l ] e ix “da 

Oir J * 


2k 


2sin(0)cos(0)— fia[n 1 B 1 e n,yi + n-,B.,e n:y ' ] e lx ' a da 
2k j 


x' x 2 '(x,y) = — — ?ia[m, A,e n '' y +m,A 1 e n ' :y +m ? A 3 e m - ,y + m 4 A 4 e nv,y ] e' x “doc 
2k i, 

+ sin(0)cos(0) — f[nfB l e n,y ' + n 2 B-,e n ’ V| ] e' x,a da 
2k 

+ sin(0)cos(0) — fa 2 [B 1 e n,y ' +B,e n = yi ] e ixa da 
2k j ' 


— (cos(0) 2 -sin(0) 2 )— J ia[n l B,e n ' yi + n 2 B 2 e n:y ' ] e'^da 


2k 


The stresses for Material 1 are 

<(x,y) = ^J[a 2 (D 1 +yD,) + 2|a|D,]e w V xa da 
271,1 

a' v , v ) (x,y) = -— f[a 2 (D, + yD, )] e |a!y e ,x “da 
2k 

jia[|a|(D, + yD 2 ) + D 2 ] e |a|y e ,x “da 


(33) 


(34) 


NASA/CR— 1999-209166 


9 



and for Material 3 are 


(x, >0 = ^ j[a 2 (C, + yC 2 ) - 2\a\C 2 ]e^ x e h0l da 

<?*■ ' (jc, y) = — — ? [a 2 (C, + y C, )> >L V l Va (35) 

2 tt 

r^x, y) = — — ? /a[-la|(C, + yC, ) + C, ]e +i| '( ;u Va 

From (33), (34) and (35), it can be seen that there are 10 constants, D|, D 2 , C|, C 2 , 
A|, A 2 , A 2 , A4, B| and B 2 (still are functions of a in the Fourier space) which must be 
determined using 10 boundary conditions. There are eight stress and displacement 
continuity conditions: 


ct''’(,y, 0) = <rj 2 '(x,0) 
o'£{x,h) = o£'{x,h) 
t^(x,Q)=t {2 '(x,Q) 
r'; V| ( -y, h ) = r[ 2 '(x,h) 
m <m (.v,0 ) = « (2, (x,0) 
u°'{x,h) = u i2, (x,h) 
v m (Jt,0) = v l2l (.v,0) 
y iM {x,h) = v a) {x,h) 


for — 00 < x < 00 


(36) 


From (36), the other constants can be expressed in terms of B| and B 2 , which in turn 
are expressed in terms of the two unknown auxiliary functions. The remaining two 
boundary conditions come from the perturbation problem, namely, 

a (x ,0) = -p. ('at, ) for a <x x <b 

(37) 

r (|Vi (.v, ,0) = -p 2 (x , ) for a<x t <b 

Here, pi and p 2 are the traction forces on the crack surfaces. D|, D 2 , Ci and C 2 are 
found in terms of Aj, A 2 , A3, A 4 , B| and B 2 by using the stress continuity conditions of 
(36). Further, by using the displacement continuity conditions of (36), the following 
linear system can be constructed: 
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( 38 ) 



where D|, D2, Ci, C2, Cn through C44 and Ji through J 4 are functions of a derived using 
MATHEMATICA, see Shbeeb (1998). 

The system of equations (38) can be solved for Ai, (i=l,..,4) in terms of the 
unknowns F|(a) and F2(a), which are components of Jj, (i=l„4) as follows: 

A =X%A ' = l-4 (39) 

./=' (3 

Here, Q is the determinant of the 4 by 4 coefficient matrix and the Qjj are the 
corresponding 3 by 3 cofactors. 

To make use of (37), the stresses of the FGM must be formulated in (X|,yO 
coordinates as 


<;,(>, ,y, ) = ~ J or (B.e'" + B.e 1 "' ) e'"“da 


1 


+ — 
2k 



(m n sin(0) + iacos(0)) 2 A n e m,,v ] e' xa da 


(40) 


(x, , y , ) = — ■— j ia(n I B,e niyi + n 2 B 2 e n:y ' ) e ,X|,1 da 


2k 


j[^((a : + mj;)sin(0)cos(0) 

_oo n — I 


+ iam n (cos(0) 2 -sin(0) 2 ))A n e m " h ] e lxa da 


(41) 


Each term in (40) and (41) must be examined for singular behavior. Upon 
substitution of (30), (27) and (29) into (40), the first integral above can be written as 
follows: 


(o 


( 2 ) 


(x,,y,)) 


<h 



a a 


( 42 ) 


where 


1 1 
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k 


k 


j' , (x l ,t) = - for 

7t J 

— oo 

p’(x,,t) = - Jor 

71 J 


h 22 e n ' y ' -h ;i e" 2y ' 
^11^22 h 1 2 h 2 1 


ia(x,-l) 


da 


h , ,e n ' Yl — h ,-,e n - yi 
h , , h :2 -h 12 h 21 


e ,a(x ' M, da 


(44) 


Since the integrands are continuous functions of a and vanish at a=0, then any 
singularity must occur as a goes to infinity. The integrand of (43) vanishes as a— >±°°, 
while that of (44) is as follows: 


h tt e tt ' y ' - h n e’" x ' 

/ij | h~r> /?|2^2I 


* ( h " !hl ) £ ,-l“l v ' 

h\ ] ^22 ^12^21 


-> + 


Tie**' 
a 2 (K 2 + 1) 


a — > ±o° 


By subtracting and adding this asymptotic value from (44) and taking the limit as yi 
goes to zero, the following is produced: 


kp (x,,t) = — [ fa 2 ( 

77 * J 


h 1 1 h.. 


2i 


Jor ( 


h | j h 22 h ph 2 j 


h, , -h, 


a 2 (K 2 +1) 
2i 


)e" 


‘da + 


(45) 


h 1 1 h h, 2 h 2l 


a~(K, + 1) 




da + SIP] 


where 


17 

SIP = lim( f 

v i — >0 J 


2ie c 


Ja( x , — t ) 


(K 0 + 1) 


da-j 


2ie 


-ay, 


Ja(x.-i) 


J 0 (K 2 +1) 


da) 


(46) 


Let a'=-a in the first integral in (46) to obtain 


— lim fe uy ' 

(K 2 +l)y.-*o{ 


sin(a(t - x . ))da = lim ^ — - 

(K 2 + l)>,^»(t-x,)- + y- 


4 1 


(k 2 + 1) t — x , 


(47) 


and 


kp(x,,t) = — [ J a 2 ( 

n J 


h,, -h 


12 


2i 


h n h 22 -h p h 2l a : (K, +1) 


-)e 


ia(x, -t) 


da + 


ja=( 


h„ -h, 2 


h , , h 2 2 -h, 2 h 2 , a“(Kp + l) 


2i 4 

)e' a<X| ~ n da 


(48) 


(k 2 + 1) t - x 


-] 


NASA/CR— 1 999-209 1 66 


12 



Similarly, let a'=-a in (48). Note that every odd power of a becomes negative. 


kp(Xpt) = — [J{or[ 
k i 


h„-h 12 


'll “12 


+ ior[ 


H 1 1 H 22 ^12“ 21 

H„ -H 


^ 1 1 ^ 22 ^12^21 


]cos(a(t - x, )) 


11 “12 


h 12 


H M H 22 -H 12 H 2I 


^11^22 ^12^21 


( 49 ) 


4i 


a'(K, +1) 


]sin(a(t -x,))}da- 


(K, +1) t -x, 


The Hjj (i, j =1,2) have the same forms as the hjj, with the only difference being in the 
sign of the odd powers of a where they are negative. It is worth noting that: 


H 1 1 Hp f h n hp 

= conjugate ! — ^ : — : — ) 


H m H 22 -HpH 


12 21 


h j j h 27 h|,h,| 


(50) 


Repeating the same procedure for (43) by splitting the integral in the same manner as 
in (48) and substituting the result along with (49) into (42) yields the following: 


Hmlcr - (X|,y.))' = 


(i) V\ e 


jSv, 


Vt — >0 


[ f ^-dt + - [/, {t)K' n (A-, ,t)dt 

K(K 2 + \) J j-x t n J ( 


(51) 


+ - 


n 


n 

J f 2 {t)K\ 2 {x x ,t)dt 


where 


Kj' l , (x 1 ,t) = -J{a 2 [X M + conj( X ,, )] cos(a( t - x , ) ) 

0 

+ ia : [-X , , + conj( X , , )] sin( a( t — x , )) }da 

Kj 2 (x,,t) =— J{a 2 [X l2 +conj(X l2 )]cos(a(t - x,)) 

0 

-> 4i 

+ ior[-X p + conj(X p ) ]sin(a(t - x, ))|da 

ot(k 2 + 1) 


(52) 


X - ( h -- h - i ) 

" h M h 22 -h 12 h 2l 

X p = ( ~ k' 2 ) 

h l l h 22 ~ h 12 h 2i 


The same procedure can be repeated for the first part of the shear stress (41) to 
obtain 
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( 53 ) 


where 


lim(r; (x, , y, ))' = [ - - - f^^-dt + -jf, (0^ 21 (*, > Odt 

v.^o >' 2 7C{K 2 +l) J a t-X l 71 { 

1 ^ 

+ — f f-,{t)K\-,{x y ,t)dt 
n J 

a 

K", ) (x 1 ,t) = -j{ia[X 2l -conj(X 2l )]cos(a(t -x, )) 

o 

4i 


+ a[X 2l +conj(X 21 ) + - 


a(K, +1) 

K 22 (x,,t) = -J{ia[X„ -conj(X„)]cos(a(t-x,)) 
o 

+ a[X 22 + conj(X 22 )]sin(a(t -x,))}doc 
x =( n,h 22 -n 2 h 2 , 

■' h„h 22 -h 12 h 21 

X = ( n :^n~ n ]^*]2 ^ 
h , , h 22 h|2^2l 


]sin(a(t -x,))}da 


(54) 


The examination of the remaining two parts of the stress equations (40) and (4 1 ) and 
application of the asymptotic expansion and the limit as yi goes to zero yields the 
following terms: 

lim(crj 2 ’ (.v, , v, )) <2> = \f ] (t)Kl; ) (x ] ,t)dt + -[Mt)Kl; ) (x ] ,t)dt} (55) 

'i-»o 2 K J n J 


fri i h - 


lim«x^ (x„ y,))' 2 ’ =£LL[L\ f l (t)Kl; > (x ] ,t)dt + - \ fM)K\;\x„t)dt] (56) 

v,-»o 2 n J k j 


where 
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tfn’CM) = J {[K,, +con;(^i)]cos[acos(0)(r-x,)] 
o 

+ i[ — yj j + cfW 7 '(^i)]sin[o:cos( 0 )(/' - x, )] }da + 

2a n (x 2 sin(0) 2 -cos(0) 2 (/-x,) 2 ) + 4fc ll (x l (/-x,)sin(0)cos(0)) 
(a:, 2 sin(0) 2 + cos (6 fit - x, ) 2 ) 2 
2c, ,x, sin(0) + 2d u cos(6)(t - x, ) ^ 
x, 2 sin(0) 2 + cos(0) 2 (r-x,) 2 

KlW^t) = J {[Y n + «wi/(^ 2 )]cos[acos(a)(f - x,)] 

0 

+ i[-Y n + co«/(y, 3 )]sin[acos(0)(r - x, )] )da + 

2a, 2 (x 2 sin(0) 2 — cos(0) 2 (r - x, ) 2 ) + 4/? l2 (x, (t -x, )sin(0)cos(0)) + 
(x 2 sin(0) 2 + cos(0) 2 (r - x, ) 2 )' 

2c, 2 x, sin(0) + 2d n cos(0)(t — x,) 
x, 2 sin(0) 2 + cos(0 fit — x, ) 2 

Kf'ix^t) = J{[yj, +c-a«y'(^i)]cos[acos(0)(/-x l )] 

0 

+ i'[-yj , + conj(Y t , )] sin [a cos (6)(t - x, )] \da + 

2 a n (x 2 sin(0) 2 -cos(0) 2 (/-x l ) 2 ) + 4b ll (x l (r-x,)sin(0)cos(0)) 
(x 2 sin(0) 2 + cos(0) 2 (/-x,) 2 ) 2 
2c, ,x, sin(0 ) + 2d n cos (6)(t - x, ) 
x, 2 sin(0) 2 + cos(0) 2 (r-x,) 2 

Kffix^t) = j{[y; 2 +cwy'(4; 2 )]cos[o!cos(0)(/-x 1 )] 

0 

+ /[-}', +c<?/y(^ 2 )]sin[acos(0)(/ -x,)]\da + 

2 a l2 (x, 2 sin(0) 2 -cos(0) 2 (/ - x, ) 2 ) + 4/? l2 (x, (t - x, )sin(0)cos(0)) 
(x 2 sin(0) 2 + cos(0) 2 (/ - x, ff 
2 c l2 x, sin(0 ) + 2d , 2 cos (6)(t - x, ) 
x, 2 sin(0) 2 + cos(0) 2 (f-x,) 2 

where Yjj ay, by, Cjj can be found in Shbeeb (1998). 
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By substituting these expressions into (40) and (41), the final system of singular 
equation is formulated: 


K 1 + 1 

k 2 + 1 


n t K i 

(-P, (X, )) = - f ^Idt + - f /, (t)K v (JC, ,r) + - f /, (f)tf„ (JC, ,0 
n J t-x, n J ‘ ?r J 


(59) 


where. 


K 11 (x,,t) (K| l , , (x l ,t) + K}f(x 1 ,t)) 

Kp(x, ,t) = ( K {V ( x , , t ) + Kp'( X| ,t)) 

4 

K,,(x 1 ,t) = -J^(K ( ,V(x l ,t) + Ox 1 ,t)) 
4 

K :: (x l ,t) = -^^(K^(x l ,t) + K^ , (Xpt)) 


3. Solution of Singular Integral Equations 

The singular integral equations (59) with the Cauchy kernels are solved for the 
unknown auxiliary functions, fi(t) and f 2 (t), by transforming them into a system of linear 
algebraic equations. In order to obtain unique results, the following conditions need to be 
incorporated with the solution: 


jfi(t)dt = 0 


a 


i = l,2 


(60) 


The singular integral equation (59) can be solved using Gaussian quadrature. For 
example, using Lobatto-Chebyshev collocation as described in Theocaris and Ioakimidis 
(1977), we obtain the system of algebraic equations in terms of discrete unknowns g(t|<) 
in the following form: 


N n 


)w k + RJx p ) = f t (x p ) 


m * *=i h~ x v 


(61) 


j = i *=i 


where p=l,...,n, the w k are the weights, the abscissas t^ are the roots of the related 
orthogonal polynomial, and R n is the error. The abscissas are 
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(62) 


t k = cos( 


(k-\)K 
n - 1 


) 


The corresponding weights are 


K 


W, = H'„ — ■ 


2(n-l) 


;w r = 


K 


n-\ 


r - 2,...,n -1. 


(63) 


The collocation points are 


.v 


p 


- cos( 


(2/7-l)7T 

2/7-2 


p = 1,. 


(64) 


Two additional equations are needed, which are generated using (60) in the 
following form: 


H 


Co 

II 

o 

(65) 

=0 

*=1 

(66) 


By combining (61), (65) and (66), the system of equations can be represented as follows: 


[^ln.v2n \S lid — f 1 


(67) 


This system can be solved by any standard method. Formally, the unknowns are 

M=UrM (68) 

Finally, the goal is to obtain the SIF in terms of gi(t) and g 2 (t). The SIF are defined 
as follows: 


k l (a) = lim J2(a — at, )<J v v (x,,0) 

,V| 1 

/c, (b) = lim ^2(x x -b)C s ^ (x ] ,0) 
k->(a) = )r U, ,0) 

M^) = lim > /2 (a, -7>)t, |V| (.v, ,0) 


( 69 ) 
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From the principal part the expressions for gi(t) and g 2 (t), Muskhelishvili (1953), the 
following is obtained for k| (a): 


k t (a) = 


2yfl^ 


(K-, +\)yjb-d 


% g 2 (a) 


Similarly, 


k, ( b ) = — 


2 ^i Bh ,,, 

r e p g 2 (b) 


(K% + 1)V b — a 


k 2 (a) = 


2V2 H\ ft , 

' ^ gl {a) 


(K-, + 1 )Jb — a 


k,(b) = - 


^ g t (b) 


(K-, + 1 )-\fb~—~ci 


(70) 


(71) 


(72) 


(73) 


Note that a=- 1 and b= 1 when solving (68). 

The strain energy release rate (SERR) can be calculated from Erdogan and Konda 
(1994). They are 


^ 4 k{k 2 + 1) , , x2 

G, (a ) = ---- ; — k t (a) 


C, ib) = 
G, (a) — 


8)U 2 (a,0) 
n{K 2 + 1) 
&H 2 (b,0) 
x(k 2 + 1) 

8jU 2 (a,0) 


k t (b) 2 

k-,{a) 2 


(74) 


^ 7T( K-, + 1) , ,^ 2 

G-,(b) = k^(b) 

%n 2 (b,0) ■ 

where G| is the opening mode SERR and G 2 is the sliding mode SERR. The total SERR 
is expressed as 


Gjia) 


G T (b ) 


XjK^ + 1) 

8//,(a,0) 


ik^a) 2 +k 2 {a ) 2 ) 


7T(K 2 + 1) 

8/i,(/>,0) 


(k^b) 2 +k 2 (b ) 2 ) 


(75) 


The verification of the solution above is accomplished by comparing the results of 
this model with that of Erdogan and Konda (1994). In this model, h is set to °o to 
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simulate an infinite FGM plate, with various values of ye. The two models give virtually 
identical results, as seen in Table 1. For y=0 (the homogeneous case), the singular 
integral equations can be reduced to a closed form solution (see for example Tada et. al. 
(1973)) producing SIF proportional to normal and shear tractions applied on the crack 
surface. 


4. Parametric Studies 

The focus of the following parametric study is limited to investigating the influence 
of the material properties of the half planes, crack length and orientation, and thickness of 
the FGM interface on the resulting driving force as measured by the SIF and SERR. To 
accomplish this, the normalized nonhomogeneity constant yh = ln((iypi) is defined using 
(3). The range of the constant is assumed to be between -3 and 3, which includes all 
known engineering materials. A negative yh represents a problem where the bottom half 
plane is stiffer than the upper half plane. A positive yh represents a problem where the 
upper half plane is stiffer than the lower half plane. Hence, if the shear moduli of all 
three phases are normalized with respect to pi, the equivalent variation of the shear 
modulus of the upper half plane takes values of approximately 0.05 to 20 times the lower 
half plane shear modulus. 

In this study, all the cases were considered under plane stress conditions with 
Poisson’s ratio = 0.3 and the materials were subjected to far field normal stress in the y 
direction. The length of the crack is chosen to be 2c, and the thickness of the interface is 
h. All geometrical dimensions are normalized with respect to c or h. Results are 
presented for the normalized mode-I and mode-II SIF, i.e., k]/ko and ko/ko, and 
normalized SERR, i.e., G|/G 0 and G 2 /Go, where k 0 = a yy (c) I/2 and G 0 = ko 2 rt(K+l )/ 8pi. 

In the first study, consider the influence of the thickness of the interface h/c and non- 
homogeneity constant yh = ln(p 2 /pi) on a crack inclined at 30 degrees such that the center 
of the crack is always kept in the middle of the interface. The distances a and b from the 
crack tips are the same from the lower and upper half plane, respectively. 

Figures 3 and 4 show mode I and II of the normalized SIF at crack tip a versus the 
non-homogeneity constant ln(pypi). Observe that as In(pypi) increases, both k ( and k 2 
decrease. The strongest effect is observed for the smallest thickness of the interface, 
while h/c = 100 may be considered as an infinite FGM plate, for which the SIF are 
virtually constant. When ln(pVfli) = 0, the plate is homogeneous, so the influence of the 
thickness of the interface disappears and the SIF become the same as for the infinite 
FGM plate. 

Figures 5 and 6 represent mode I and II normalized SIF at crack tip b versus 
ln(p 3 /pi). The SIF curves increase with increasing lnlp^/pi), which is different from the 
behavior at tip a, except for the case of the infinite FGM plate. In addition to this 
behavior, the magnitude of the SIF tends to be higher at crack tip b than at tip a, 
especially for extreme values of ln(|i 3/(11). 

Modes I and II of SERR are shown in Figure 7 and 8 for crack tip a and in Figure 9 
and 10 for crack tip b. Notice that at both tips, SERR are decreasing with increasing 
stiffness of the upper half plane. Specifically, SERR at crack tip b behaves differently 
from the SIF at this tip. The behavior of SERR is more physically intuitive than the 
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unexpected behavior of the SIF. It should be recalled that the SERR is calculated using 
SIF and local material properties, so it contains more information than the SIF. For this 
reason the remaining parametric studies are discussed using only SERR data. 

The influence of the orientation angle 0 on the relation between SERR and the non- 
homogeneity constant at crack tip a for the case of the interface thickness h/c=2 under 
uniform shear stress at infinity is shown in Figures 1 1 and 12. As expected, the highest 
mode-I SERR is obtained for the smallest angle because of the high normal traction 
component acting on the surface of the crack. It can also be expected that the highest 
magnitude mode-II SERR is produced for 0=45 degrees, since the shear traction 
component is maximized then. The behavior of the SERR at crack tip b is similar to the 
behavior at crack tip a, as can be deduced by comparison of Figures 7 and 9, and of 8 and 
10. Hence, the magnitudes of the mode-I SERR at crack tip b are higher than at crack tip 
a, while the magnitudes of the mode-II SERR are smaller at crack tip b. In fact, the total 
SERR at crack tip b is equal to the total SERR at crack tip a. 

In the next study, assume that the crack orientation is 30 degrees from the horizontal, 
and the thickness of the interface is h= 1 . Crack tip a is fixed at the distance a/h = 0. 1 
from the origin while crack tip b is at distance equal to b/h= 0.3, 0.7 and 1.1 along the xi 
axis, making the half of the crack length c/h = 0.1, 0.3 and 0.5, respectively. Modes I and 
II normalized SERR at crack tip b are shown in Figures 13 and 14, respectively. Notice 
that the crack length significantly changes both SERR modes in the case of negative non- 
homogeneity constant. For the case where the upper half plane is stiffer than the lower 
half plane, the longer crack produces smaller normalized SERR. 

Finally, assume constant crack length, constant orientation at 30 degrees and 
constant thickness of the interface FGM h/c = 2, and examine the influence of the 
position of the crack along the Xi axis. Figures 15 and 16 show modes I and II normalized 
SERR at crack tip a versus the non-homogeneity constant for the crack defined by tip 
positions varying from 0.2 to 1.6 from the origin. Notice that the largest modes I and II 
SERR are obtained when a/c = 1.6 for negative non-homogeneity constant and when a/c 
= 0.2 for positive lndii/pi). Hence, the closer the crack tip is embedded to the stiffer 
material, the smaller the normalized SERR. 


5. Conclusions 

The analysis of an arbitrarily oriented crack in a strip of FGM sandwiched between 
two isotropic homogeneous half planes is done using singular integral equations. The 
equations are solved using Lobatto-Chebyshev integration, and give accurate results for 
mix-mode SIF and SERR. 

Parametric studies show that SERR contain relevant information that is missing in 
the SIF, and therefore it is recommended that SERR be used as a driving force parameter 
for fracture problems of a crack in FGM. The model has shown that SERR are sensitive 
to the ratio of the shear moduli used as non-homogeneity constant. They are also 
sensitive to the ratio of thickness to the crack length. The longer the crack or the thinner 
the interface, the larger the SERR produced for negative values of ln((iV(ii), and the 
smaller the SERR produced for positive values of ln(|i ? /(i|). 
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The crack orientation influence shows that SERR is proportional to the traction 
forces at the crack surface assumed in the perturbation problem. As the lower half plane 
becomes stiffer (ln((j.3/(ii) becomes more negative), both modes of the SERR become 
larger for every crack orientation. Clearly, the proper selection of the FGM parameters 
can reduce the driving forces of a crack embedded in the interface material. 
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Table I Verification of the Solution. 


yc 

Konda and 
Erdogan (1994) 
k,(a)/Vc 

Present Study 
k a (a)/ Vc 

Konda and 
Erdogan (1994) 
k 2 (a)/ Vc 

Present 
Study 
k2(a)/ Vc 

0.25 

1.036 

1.036 

0.065 

0.062 

0.50 

1.101 

1.101 

0.129 

0.122 

1.0 

1.258 

1.260 

0.263 

0.243 


c=(b-a)/2 
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APPENDIX A 


EXPRESSIONS OF THE CONSTANTS 


It should be pointed out that all the algebraic manipulation were either verified 
done by MATHEMAT1CA®. 

D| = A, + A 2 + A 3 + A 4 -4- g j Bj + g 2 B 2 

D, = (m, — jot|)A, + (m, -|oc|)A, + (m, -|a|)A, +(m 4 — |cx|)A 4 

-(|a|g, - — )B, -(|a|g 2 )B 2 

11 a a 

C, =e |ot|h [(l-h(m 1 +|a|))A,e m ' h +(1 -h(m 2 +|a|))A 2 e m2h 
+ (1 - h(m, +|a|))A,e" l,h + (1 - h(m 4 + |a|))A 4 e"’ jh 

+ e i«hTan,e) { ( (1 _ h|a|)-^^-)B'e cn ' hscc(ei + (g, (1 - hlal) -^^-)B'e cn:hsec<e> } ] 

a a 

C 2 = e |a|h [(m, + |oc|) A ,e m,h + (m 2 + |a|)A 2 e msh + (m, + |a|)A 4 e ni ’ h 
+ (m 4 +|a|)A 4 e mjh +e iahTan(e, {(|a|g, + ^)B;e cn ' hsec(e ' 

+ (lalg , + ^-)B'e cn2hsec<0) }] 
a 


cos(0) 
gl ~ ” 


(cn, — iacos(0) 2 ) 2 


or 

COS(0), 

g, = — (cn, -iacos(0) ) 

O' 

g, =cos(0) 2 (asin(0)-icn l )(acos(0) + icn, sin(0)) 
g 4 = cos(0) 2 (asin(0) — icn,)(acos(0) + icn 2 sin(0)) 
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C n =2^m,(ic 1 + l) + a(v, -2v, (K ~- + 1) 

a ‘ a 


C p =2 — /Mk, + l) + a(K:, — 2 k-j _ i>_ 

a ' " a 

led /?Z?(JO, + l) 

C n = 2 — Ytt^ ( k*j + 1) + Of( /c 9 — 2/c, — 1) 1 — 

a ‘ a 

C )4 = 2— /m 4 (k-, + l) + a(K\ -2k ] - \)- ^—2 + ^ 
a ‘ a 


C 2] = 2m, (K t — 1) — 2|or|(Kr, + 1)- 
C 22 = 2m 2 (K ] - l)-2\a\(K t +1)- 


(k 2 - 3 )m, 2 -orCic, + 1 ) 
m, -7 

( k :, -3 )m; -a 2 (K^ + 1 ) 


m 2 - 7 

l "> . 


C 2 , = 2m, (k*, — 1> — 2|of|(KT, + 1)- 


C 24 = 2 m 4 (K’ 1 - l)- 2 a(jc l + 1 ) 


(k 2 - 3 )m; -ot(k 2 + 1) 

n h - y 

( k 2 -3 )m 4 -cr(K\ + 1 ) 


m 4 - 7 


a 


C 31 = [(2a(l-/c ? )-2— m,(l + K‘ 4 )) + a(/c, -3)- 

a 

lal 


C,, =[(2ff(l-if,)-2Q m, ( 1 + v , )) + of(y, - 3) - v " - ■ - — ]e 

a ‘ 


a 


C„ =[(2a(l-je,)-2 — m,(l + K\)) + a(K\ - 3) - 

a 

C, 4 = [(2a(l-K\,)-2 — m 4 (l + x\)) + a(K\ -3)- 

a 


m 2 (K 2 

+ 1) 

a 


m\(K 2 

+ 1) 

a 


m 2 (K 2 

+ 1) 

a 


m 2 ( k 2 

+ 1) 


a 


}e 


]e 


(k, -l)m 2 -cr(K% + 1) b 
]e 


C 41 = [(4|a| + 2(^3 - l)(m, + |a| 
C 42 = [(4|a| + 2(/c, - l)(m, + |a| 


/n, - 7 


( k 2 - l)m; -«-(k- 2 + 1 ) 1 
m 2 - 7 


-k" 


C 4 , = [(4|a| + 2<ir, - l)(m, + |a|» - ^^k"'* 


c« = [(4|a| + 2(k, - 1 )(m 4 + |a|)) - 


m, - 7 
m 4 — a 2 

'»4 - Y 


( y,-l)m 4 --or(y 2 + l) ^ 


»f]/i 


w->// 


hi \h 


m±h 
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n 

7 , = 5 , | /, (t^'e^^dt + R x J f 2 (t)e p, e~ iCttcos(e) dt 


S, = 


g 6 i7? 2l - gjChyj . 2a (1 - g, )(g 2 c/; 2 , - g,c/i 22 ) | 2/(ir, + 1 )(g 4 cft 2 , - g 3 cft 22 ) 

|a|(c// u c/j 22 - ch n ch 2 , ) 


ch u ch 

22 - ch x2 ch 2x 


ch xx ch 22 -ch x2 ch 2x 


R , = 


, 2«(l->c,)(g|CA, 2 -g 2 c/t, l ) | 2i'(KT| + l)(g 3 c/t, 2 -g 4 c/»,,) 

|a|(c/i n c/i 22 - ch l2 ch 2x ) 


ch xx ch 22 - ch x2 ch 2x 


ch xx ch 22 - ch xl ch 2x 


p n 

J 2 = S.jUt^'e-^dt + R 2 J f 2 (t)e p 'e~ MBl dt 


S Si ch 22 ~ 8 »E!hl I 2 H^ + ’Si ~ ^ , 2t '( y i - 0 (g 4 c ^ 2 l ~ g ^22 ) 

ch xx ch 22 —ch x2 ch 2x ch u ch 22 -ch x2 ch 2x a(ch xx ch 22 — ch x2 ch 2x ) 

D ^ 8 c7j , 2|a|( 1 + v, X^c/?, , -g,c7i 12 ) 2i(ir l - l)(g,c/? P -g 4 ch xx ) 

= — : : 1 h 


ch xx ch 22 -ch x2 ch 2x 


ch xx ch 22 — ch r ch 2x 


a(ch xx ch 22 - ch n ch 2 , ) 


It l> 
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APPENDIX B 
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Figure 1. Formulation of the Perturbation Problem. 

(a) . The Original Problem. 

(b) . The Elasticity Problem. 

(c).The Mixed Boundary Value Problem. 
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(b) (c) 

Figure 2. Methodology of the solution of the Perturbation Problem 
(a). The Mixed Boundary Value Problem 
(b). Infinite FGM Strip Without Crack 
(c). Infinite FGM Plate With Crack 
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Figure 4. Normalized mode II SIF at crack tip (a) for various h/c, 0=30 deg. and 
center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 5. Normalized mode I SIF at crack tip (b) for various h/c, 9=30 deg. and 
center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 6. Normalized mode II SIF at crack tip (b) for various h/c, 0=30 deg. and 
center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 7. Normalized mode I SERR at crack tip (a) for various h/c, 0=30 deg. and 
center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 8. Normalized mode II SERR at crack tip (a) for various h/c, 0=30 deg. 
and center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 10. Normalized mode II SERR at crack tip (b) for various h/c, 0=30 deg. 
and center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 1 1 . Normalized mode I SERR at crack tip (a) for various 0, h/c=2 and 
center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 12. Normalized mode II SERR at crack tip (a) for various 0, h/c=2 and 
center of the crack is located at h/2, under loading of uniform normal stress. 
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Figure 13. Normalized mode I SERR at crack tip (b) for 0=30 deg., same h and 
fixed crack tip (a) and movement of crack tip (b), under loading of uniform 
normal stress. 
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Figure 14. Normalized mode II SERR at crack tip (b) for 0=30 deg., same h and 
fixed crack tip (a) and movement of crack tip (b), under loading of uniform 
normal stress. 
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Figure 15. Normalized mode I SERR at crack tip (a) for 0=30 deg., h/c=2, 
constant crack length and various positions of crack, under loading of uniform 
normal stress. 
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Figure 16. Normalized mode II SERR at crack tip (a) for 0=30 deg., h/c=2, 
constant crack length and various positions of crack, under loading of uniform 
normal stress. 
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